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Abstract We study superfiuid behavior of a gas of spatially indirect magnetoex- 
citons with reference to a system of two graphene layers embedded in a multilayer 
dielectric structure. The system is considered as an alternative of a double quantum 
well in a GaAs heterostructure. We determine a range of parameters (interlayer dis- 
tance, dielectric constant, magnetic field and gate voltage) where magnetoexciton 
superfiuidity can be achieved. Temperature of superfiuid transition is computed. 
A reduction of critical parameters caused by impurities is evaluated and critical 
impurity concentration is determined. 
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1 Introduction 

Recent progress in creation of heterostuctures with two graphene layers separated 
by a thin dielectrics [1 opens possibilities to use graphene for creation of multiple 
quantum well structures with separately accessed conducting layers. In [1] Si02 
substrate and AI2O3 internal dielectric layer were used. Another promising dielec- 
tric is hexagonal BN It has a number of advantages, such as an atomically 
smooth surface that is free of dangling bonds and charge traps, a lattice constant 
similar to that of graphite, and a large electronic bandgap. 

The attention to graphene heterostructures is caused, in some part, by the idea 
to use them for a realization of superfluidity of spatially indirect excitons [3ll31[5l 
|6j[71[8l|9] . Bound electron-hole pairs cannot carry electrical charge, but in bilayers 
they can provide a flow of oppositely directed electrical currents. Therefore, exciton 
superfiuidity in bilayers should manifest itself as a special kind of superconductivity 
- the counterflow one, that means infinite conductance under a flow of equal in 
modulus and oppositely directed currents in the layers. 
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The idea on counterflow superconductivity with reference to electron-hole bi- 
layers was put forward in [lOllllj . The attempts to observe counterflow conduc- 
tivity directly were done [12lll3|[Ti] for bilayer quantum Hall systems realized in 
GaAs heterostructures. In the latter systems superconducting behavior might be 
accounted for magnetoexcitons [ISIfTB] . The effect is expected for the filling factors 
of Landau levels Vi = 2Tv£^ni {£ = ^Hc/eB is magnetic length, is the electron 
density in the i-th layer) satisfying the condition 1/1 + 1/2 = 1- The role of holes is 
played by empty states in zero Landau level. In experiments |12lll3l[Ti] an expo- 
nential increase of the counterflow conductivity under lowering of temperature was 
observed, but zero-resistance state was not achieved. The latter can be explained 
by the presence of unbound vortices [17, 18, 19 . Such vortices may appear due to 
spatial variation of the electron density caused by disorder. 

To demonstrate counterflow superconductivity quantum Hall bilayers should 
have the parameters that satisfy two additional conditions: d < £ and ^ < a^, 
where d is the interlayer distance, and a*g = ef? je^vn* is the effective Bohr radius 
(e is the dielectric constant of the matrix, and m* is the effective electron mass). 
The first inequality comes from the dynamical stability condition. For balanced 
bilayers (v\ = 1/2) the mean-fields theory yields d < 1.175^. The second inequality 
is the condition for the Coulomb energy e^/e£ be smaller than the energy distance 
between Landau levels. In GaAs a*g ^ 10 nm and the condition ^ ^ a|j is fulfilled at 
rather strong magnetic fields B > 6 T (actually, the experiments |121ll3|[Ti] were 
done at smaller fields). At d < 10 nm the interlayer tunneling is not negligible 
small and may result in a locking of the bilayer for the counterflow transport at 
small input current [201121] . At larger input current the system unlocks, but the 
state becomes nonstationary one |221l23ll2i] that is accompanied by a dissipation 
(the power of losses is proportional to the square of the amplitude of the interlayer 
tunneling p^ll^g] ). 

The idea to use graphene for the realization of electron-hole superfluidity in 
quantum Hall bilayers [6l|71[8l|9] looks very attractive. The distance between Lan- 
dau Levels in monolayer graphene is proportional to the inverse magnetic length, 
magnetic field does not enter into the condition of smallness of the Coulomb en- 
ergy, and small magnetic fields can be used. Smaller magnetic fields correspond 
to smaller critical temperature, but, at the same time, they correspond to larger 
critical d. Use of large d allows to suppress completely negative effects caused by 
interlayer tunneling. 

In this paper we concentrated on three questions. First, we determine, in what 
range of internal parameters and external fields magnetoexciton superfluidity can 
be realized. Second, we evaluate critical temperature for pure system. Third, we 
consider its reduction caused by electron-impurity interaction. Our study extends 
the results of [8^, where a system of two graphene layers embedded into a bulk 
dielectric matrix was considered. Here we investigate structures with one and two 
graphene layers situated at the surface. 



2 Conditions for the electron-hole pairing in zero Landau level 

Quantum Hall effect in graphene is characterized by unusual systematics of Landau 
levels and the additional four-fold degeneracy connected with two valleys and 
two spin projections [25]. The energies of Landau levels in graphene are E±pf = 
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Si^^2\N\, where N = 0,1,2,..., and vp ^ 10^ m/s is the Fermi velocity. In a 
free standing graphene the A*" = Landau level is half-filled. A state with only 
completely filled Landau levels corresponds to a plateau at the Hall conductivity 
plot (dependence of axy on electron density). A free standing graphene is just 
between two plateaus [2^ . A given quantum states in zero Landau is characterized 
by the guiding center index X and the combination of the spin and valley indexes. 
Below we call four possible combinations, the components, and numerate them by 
the index l3 = 1,2,3,4. 

We describe electron-hole pairing in zero Landau level in graphene by the wave 
function that is a generalization of the wave function [151 to the multicomponent 
case 

f) X 

Here c^^x electron creation operator (the operator that fills a given state in 

A'' = Landau Level), |0) is the state with empty zero level, i is the layer index. 
The u — V coefficients satisfy the condition jtt^p + =0. The function ([1]) can 
be rewritten in the form 

i'^) = nn("'3 +«/3c+ (2) 

where h'^i^x ~ '^ifix is the hole creation operator, and the vacuum state is defined 
as \vacl) = Yip Ylx ^1/3x1'-') • '^^^ that the function ([2]) is an analog of the 
BCS function in the Bardin-Cooper-Schrieffer theory of superconductivity. 

The quantity up = |tt^P — k/9p gives the filling factor imbalance for the com- 
ponent /?. The order parameter of the electron-hole pairing reads as Zi^ = u^vp = 



1 — u'^e^^ /2. If a given component is maximally imbalanced {up = ±1) the order 
parameter is equal to zero. 

If a one component bilayer system is balanced, the order parameter for the 
electron-hole pairing is maximum. But if the number of components is even, the 
balance up = Q can be reached &t up = \ for half of the components and 
Up = —1 for the other half. In the latter case all = 0. As is shown below, 
just such a state corresponds to the energy minimum. In other words, in balanced 
graphene bilayers electron-hole pairing does not occur. 

At nonzero imbalance '^pUp ^ 0,±2,±4 at least for one component up ^ 
±1, and electron-hole pairing may occur. Nonzero imbalance can be provided by 
electrical field directed perpendicular to the layers. Such a field can be created by 
a voltage difference applied between top and bottom gates (see fig. [T|) . 

We consider the general structure "dielectric 1 - graphene 1 - dielectric 2 - 
graphene 2 - dielectric 3" with three different dielectric constants ei, £2 and £3. 
Dielectrics 1 and 3 are assumed to be thick (much thicker than the distance be- 
tween graphene layers d). Solving the standard electrostatic problem we obtain 
the Fourier components of the Coulomb interaction Vnr for the electrons located 
in i and i' graphene layers 

V ( ) = ^I^ £2 +£3 + (£2 -£3)e~^'^'' , . 

"^'^^ q (£2+£3)(£2+£l)-(£2-£3)(£2-£l)e-2f^' 
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Fig. 1 Schematic view of the system under study. CI - C4 are the contacts. 



V22(q) 



47re 



£2 +£l + (£2 -£l)e 



-2qd 



(£2 + £3)(£2 + £l) - (£2 - £3)(£2 - £l)e-29'' ' 



£26 



q (£2 + £3)(£2 + £l) " (£2 - £3)(£2 - £l)e ' 



(4) 
(5) 



For electrons in TV = Landau level in graphene the Hamiltonian of Coulomb 
interaction has the form 



X,X' P,/3' q 

'4l3X+qyt^ /2^t' 13' X' -QyP /2'^l' X' +qyP /2C-lpX-qyP 



(6) 



where S is the area of the system. The interaction with the gate field is described 
by the Hamiltonian 



^ E (^tajfCi/SJC - C+3xC2/3Jf ) , 



(7) 



xp 



where Vg is the interlayer voltage created by the external gate (bare voltage). 
Rewriting the wave function ^ in the form 

1"^) = n n ('^"^ T^^"^^ + t4^x) io>, (8) 

X fj ^ ' 

and computing the energy in the state (|8]) we obtain 
S 



mf 



8tvP 



W E COS cos 0(3' -JoJ2 ~ (2e Vg + ) E cos 6^^ J , (9) 



where W = e^d/e2^^ is the energy of direct Coulomb interaction. The exchange 
interaction energies 

determine the parameters Jq = ( Jii + J22)/2 — J12 and Jz = Jii — J22- The relation 
between 9p and i/^ is given by equation = coaOp. 

Taking into account the inequalities W > Jo, and Jn, J22 > J^i2 (that can 
be checked directly) we find that at Vg = the minimum of (|9]) is reached at 
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i/1 = i/2 = 1, i's = i'4 = — 1- It indicates the absence of electron-hole pairing in 
balanced systems. 

If Vg ^ and belongs to one of the intervals 

nW + J22 - J12 < eVg < (n + 2)W - Jn + J12, (10) 

where n = —4, —2, 0, 2, the energy minimum is reached at i/^^ ^ ±1 for one of the 
components. We will call such a component the active one. 

Let us, for instance, consider the interval (|10p with n = 0. Then the energy 
minimum is reached at 

eVg + Jz/2-W 



W-Jo 



The case = (with maximum order parameter) corresponds to the voltage 

eVg = ^^ + W. (11) 

Eq. (|lip determines the relation between magnetic field and the gate voltage Vg. 
To keep up^ = the gate voltage should be varied synchronically with B. In 
particular, at Jz = (ei = £3) the quantities Vg and B are linearly related: 

Vg = (12) 

£2 

where a ^ 1/137 is the fine structure constant (the relation (|12p is given in SI 
units) . 

If only the gate voltage or magnetic field is varied, the order parameter (and 
the critical temperature) changes nonmonotonically reaching the maximum at the 
point determined by pip . 



3 Collective mode spectrum and phase diagram 

The components that belong completely to one layer do not take part in the 
pairing. In what follows we consider the dynamics of only the active component. 
We describe the active component by the wave function 

1"^) = n (^^^ 'f<X+QyP/2 + e'(«^^+'^-) sin ifcl^_^^,,A 10) (13) 
X ^ ^ 

(here and below we omit the component index). Eq. (|13p describes the state with 
nonzero counterfiow currents. To illustrate this statement we neglect for a moment 
the order parameter fluctuations (i^x = 0, = ^a). 
The order parameter is determined by the equation 

A(v) = J2 rx{r)i'X'{r){^\c+xC2.x'\'P)- (14) 
X,X' 

where 

, , 1 (^-x)^ 
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is the single-particle wave function in the coordinate representation, Ly is the 
width of the system. Substitution (|13|1 into (|14p yields 

A(r) = — 2— e ^ . (15) 

One can see from Eq. (|15p that Q = {Qx,Qy) is the gradient of the phase of the 
order parameter. 

Computing the energy in the state (|13|) and neglecting the fluctuations we 
obtain 

Eo = ^ ([W - Fs{Q)] cos" 9a - Fd{Q) sin" 9a) , (16) 

where 

Fs{q) = ^J^ pMpqe^)[Vii{p) + V22ip)]e- — dp, (17) 

and 

Foiq) = TT / pMpqe^)Vi2{p)e-^dp. (18) 

Electrical currents can be found from a variation of the energy caused by a 
variation of the vector-potential 

5E = ~-^j cf^Y,i^5J^^■ (19) 

Here is the in-plane component of the vector-potential in the layer i. To obtain 
the explicit expression for the variation p9p we replace the phase gradient in 
p6p with the gauge-invariant expression Q — -^{Api i — Ap; 2)1 where A^; j is the 
parallel to the graphene layers component of the vector potential in the layer i. 
Then, using (|19|) one finds the currents 

e sin" 9a dFpjQ) 



At small gradients Q£ <^ 1 Eq. (|20|) is reduced to 

g 

ji = j^PsoQ, (21) 
where coefficient of proportionality between the current and the phase gradient 

/)2 fOO 2,2 



Pso = iTT^sin^ Oa p^Vi2{p)e-— dp (22) 

is called the zero temperature superfluid stiffness (the definition is given in the 
next section). Since we neglect fiuctuations, the expression (|20p yields the current 
at T = 0. 

Implying the fiuctuations of the amplitude and the phase of the order param- 
eter are small one can present the energy as 

E = Eo + E2 + ... (23) 
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The quadratic in fluctuations term can be diagonalized: 

E2 = ^[m;,{~q)Kzz{q)mz{q) + i(/3(-g)-ft:^^(g)<^(q) 
Q 

--^{imz{-q)Kz^{q)ip(q) + c.c.)], 

where 



(24) 



M) 



1 /27rP 



S 



(cos 6x — cos 6a) e 



-iqX 



X 



2nP 



-iqX 



(25) 



X 



are the Fourier components of the fluctuations. 

Eqs. (|24p yields the energy of fluctuations with the wave vector directed along 
the X axis. The component of the matrix K can be presented in form independent 
of the choice of the direction of the coordinate axes 



where 



Kzz(,Ci, Q) = i/(q, Q) - Fs(iql) + F,3(|Q|) + H(q, Q) COt" Oa, 

^^^^(q.Q) =sin^^aH(q,Q), 
KzM, Q) = - cosea[FD{\c\ + Ql) - Fo(|q - Q|)]/2, 

"V^ll(q) + T/22(q) 



(26) 



^^(q,Q) 



S(q,Q) 



Vi2(q)cos {\ci X Cl\f 

fg(|q + QI) + fg(|q-QI) 



(27) 



The quantities K^jj (g) in (|24p are expressed in terms of (|26|) as AT^^ (g) = JCq,^ (q, Q) 

The quantity h cos 9x /2 can be treated as a ^-component of the pseudospin and 
it is canonically conjugated with the phase i^x- The Fourier transformed quantities 
(|25p are deflned as canonical variables as well. The equations of motion for the 
quantities mz{q) and ^{q) read as 



h 



dt 



h 



dmzjq) 



dt 



2Kzz{q)mz{q) - iKz^{q)f{q), 
7;Kvv{q)ip{q) - iKz^p{q)mz{q). 



(28) 



Eqs. (|28p yield the collective mode spectrum i?(g, Q) = \J Kipip{q)Kzz{q) +Kzip{q). 
Rotating the axes one obtains the excitation spectrum at general q 



^2(q,Q) = VK^tl, Q)Kzz{ci, Q) + Kz^ici, Q). 



(29) 



At Q = the spectrum (|29|) is isotropic. It can be presented in the Bogolyubov 
form 

(30) 



In Eq. ([30] 



£g = Fd{0) - Foiq) 



(31) 
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Fig. 2 Graphene hererostructures under study. 

is the kinetic energy (e^ ^ f?q^ jlM at qi ^ 1, where M is the magnetoexciton 
mass, see, for instance [27 ), and 

7, = [H(q,0) - + Fo{q)\ Sin" So (32) 

has the sense of the exciton-exciton interaction energy (that includes the direct 
and exchange parts). 

The condition for the dynamical stability of the state (|13p is the real valueness 
of the excitation spectrum (|30p . This condition determines the diapason of djl and 
Ei where superfluid magnetoexciton state can be realized. To be more concrete we 
consider three types of heterostructures. Type A is a graphene-dielectric-graphene 
sandwich with two graphene layers at the surface, type B is a graphene-dielectric- 
graphene-dielectric structure with one such a layer, and type C is a system of two 
graphene layers embedded in a dielectric matrix (Fig. [2]). For simplicity, we imply 
the same dielectric constants e for the interfacial layer and the substrate. 

The dynamical stability condition is fulfilled at < dji < rfc(e), where dc(e) 
depends on the imbalance parameter vp^ = Va- The dependence dc{e) aX Da = Q is 
shown in Fig. |3l 

The requirement for the Coulomb energy be smaller than the distance between 
Landau levels yields the restriction on e. Since we study the pairing in TV = 
Landau level we compare the Coulomb energy with the energy distance between 
TV = and TV = 1 levels ojc = ^Plhvpli. 

We have four parameters that characterize the Coulomb energy VK, Jn, J22 and 
J12. At dji < dc the largest of them is Jn (the intralayer exchange interaction in 
the graphene layer at the surface) . Therefore, it is natural to consider the condition 

Jn < ujc (33) 

as the additional restriction on the parameters. Eq. (|33p can be rewritten as e > 
ec{d/l). The quantity Ec can be understood as a critical dielectric constant. The 
dependence ec{d/l) is also shown in Fig. [S] 

Two conditions d/(. < dc{e) and e > ec{d/£) determine the range of parameters 
where one can expect a realization of electron-hole pairing and magnetoexciton 
superfluiduty in graphene bilayer systems. 
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4 Critical temperature 

In a bilayer graphene heterostructure with a fixed d the magnetoexciton superflu- 
idity can be realized in a wide range of magnetic field. Variation of B at fixed gate 
voltage results in a change of imbalance of the active component. Simultaneous 
tuning of Vg allows to keep zero imbalance Va = and maximum order parameter 
under variation of _B. In this section we study the dependence of critical temper- 
ature on magnetic field implying such a simultaneous tuning. 

Superfiuid transition temperature is given by the Berezinskii-Kostelitz-Thouless 
equation [15] 

Tc = ^Ps{Tc), (34) 

where ps (T) is the superfiuid stiffness at finite temperature. The superfiuid stiffness 
is defined as the coefficient in the expansion of the free energy in the phase gradient 
F = Fo+f d^rps(V((9)^/2. In a weakly nonideal Bose gas it is equal to ps = fi^ns/m, 
where Us is the superfiuid density. As was shown in previous section, superfiuid 
stiffness determines also the supercurrent. 

Taking into account linear excitations we present the free energy F = Eq — TS 
in the following form 



In h - e ~] 

q 



(35) 
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Fig. 4 Critical temperature vs magnetic field for A, B, and C structures. Temperature is given 
in units of = jsd, magnetic field, in units of B^i = </</7rd^. 



Expansion of Eq. (|35p yields the following expression for the superfluid stiffness 



Ps 



(T) = p,o + ^ ^ 



dQ2 



Na 



Q=0 



r 



Q=0 

(36) 

It follows from (|36p and (|29p that Ps{T) < pso (thermal fluctuations reduce the 
superfluid stiffness). 

For the spectrum fi(q) = E{q) + hqy (where v = h\7(p/m is the superfluid 
velocity) (|36|) yields the well-known answer for the superfluid density |28]. Eq. 
(|36|) generalizes the results [28| for the general case. 

The dependence of critical temperature on magnetic field at Da = and e = 4 
is shown in Fig. 3] One can see that the maximum critical temperature is reached 
approximately at B ~ 0.5B^, where = 4>/ttcP with 4> = hc/2e, the magnetic fiux 
quantum. 



5 Influence of impurities on the critical parameters 

In the previous section we have determined the infiuence of thermal fiuctuations 
on the superfluid stiffness. In this section we consider the effect of reduction of the 
superfiuid stiffness caused by the interaction of magnetoexcitons with impurities. 

The Hamiltonian of the interaction of the active component with impurities 
can be presented in the form 

H^mp = ^J2^z{c^){Mn)-p2{c^)), (37) 

q 

where Uz{ci) = t^i(q) — t^2(q), t^i(q) is the Fourier-component of the impurity 
potential in the layer i, and 

ft(q) = J2 ""tx+iM^^x-i^ '^''P [-iqxX - — j (38) 
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is the Fourier component of the electron density operator for the active component. 

In the state (|13p the energy of interaction with the impurities expressed in 
terms of mz{q) reads as 

Ejmp = const + Uz{q)mz{q), (39) 
q 

where 

The interaction (|39p induces the fluctuations of the density and the phase of the 
order parameter. Their values can be obtained from the Euler-Lagrange equations 



Smz{q) 
6Lp{q) 

where E is the energy of the system, described by the Hamiltonian H = He + 
Hq + Himp in the state (fT3)) . 

Eqs. (|40p solved in linear in impurity potential approximation yield 

'"'^''>- 2Kzz{q)KMl)'KU<iy ^ ' 

^^"^ Kzz{q)K^A<l) - Kl^ilY ^ 

Substituting (|4ip . (|42p into the expression for the energy one finds the correction 
to the energy caused by the electron-impurity interaction 

^^ = -7yf^4'?)f^-(-9)T^^T#^rx^^rTT- (43) 

4Y Kzz{q)K^Aq)~ K^zv{q) 

In Eq. (|43p the contribution of fluctuations with the wave vectors directed along x 
is taken into account. Summing the contribution for all wave vectors one obtains 

For simplicity, we specify the case where impurities are located in graphene 
layers. Then the Fourier-component of the impurity potential can be presented in 
the form 

Uz{^) = Y.e''"'^Uz,M), (45) 

a 

where Va are the impurity coordinates, and iiz,i(q) = 'Ui,i(q) — U2,i{f\) with Mfe_i(q), 
the potential in the layer fc of a single impurity centered at r = in the layer i. 
Averaging over impurities yields 

le 2 



(46) 
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Fig. 5 Critical impurity concentration vs magnetic field for charged impurities located in 
graphene layers. 



where rij^p is the impurity concentration in a layer. 

At Q£ ^ 1 the energy (|46p can be expanded in series as 



AE = AEo + ^ApsQ^, 



(47) 



where 



"'imp 
8nPS 



E 



(Kl(q)|^ + K2(q)|^)e-^ 



I 



aQ2 I 



^ I — m — 



V 



(48) 



is the correction of the superfluid stiffness. One can check that the correction Aps 
is negative. Thus, the interaction with impurities resuhs in decrease of critical 
parameters. 

At i>a = Eq. (gS} is reduced to 



_ "-imp ^ (i^z,l (q)P + |Mz,2(q)|^) e "2 



S 



i^L(q,0) 



g^Vi2(g)e ^ 
327r2 



(49) 



where pso (Eq. (122^ 1 is taken at Qa = 1^/2. 

The shift of critical temperature is evaluated as ATc/Tc ~ ^Ps/pstEl- We define 
the critical impurity concentration n^^^ as a concentration at which Aps/ pso = 1 • 
We consider charged impurities with the potential ?tz,i(q) = (— l)*(Vi2(q) — V'ii(q))- 
The dependence of critical impurity concentration on magnetic field at e = 4 and 
i>a = is shown in Fig. [S] We also evaluated critical concentrations for neutral 
impurities. These concentrations are much larger, and the influence of neutral 
impurities can be neglected. 



^ Since in our approach we assume smallness of Aps/ Pso it is just an estimate 
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6 Conclusion 

In conclusion, we present some estimates. Let us specify the type B structure (the 
one used in [Tj) with d = 20 nm and e = 4. For this structure the maximum critical 
temperature Tc « 3 K (in pure case) is reached in magnetic field B ^ 0.8 T. At 
such B the critical impurity concentration is n^^^ ^ 2- 10® cm^^. The gate voltage 
determined by Eq. is Vg ^ 6 mV, that corresponds to electrostatic field E ^ 3 
kV/cm. 

Basing on the results of our study we may state the following. 

1. Graphene bilayer structures are perspective objects for the observation of 
magnetoexciton superfluidity. The advantages are smaller magnetic flelds and no 
restriction from above on physical interlayer distance, that means the possibility 
to suppress completely interlayer tunneling. 

2. Gate voltage should be created between graphene layers for a realization of 
magnetoexciton superfluidity. 

3. Certain conditions on dielectric constant and on the ratio between interlayer 
distance and magnetic length should be satisfied. 

4. Structures with graphene layers situated at the surface have larger critical 
parameters. 

5. Neutral impurities are not dangerous for the magnetoexciton superfluidity, 
but the concentration of charged impurities should be controlled. 

Acknowledgements This study is supported by the Ukraine State Program " Nanotechnolo- 
gies and nanomaterials" Project No 1.1.5.21. 
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